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Abstract. Bright plane soliton solutions of an integrable (2+1) dimensional (n + 1)- 
wave system are obtained by applying Hirota's bilinearization method. First, the 
soliton solutions of a 3-wave system consisting of two short wave components and 
one long wave component are found and then the results are generalized to the 
corresponding integrable (n+ l)-wave system with n short waves and single long wave. 
It is shown that the solitons in the short wave components (say 5^^^ and S'^^^) can 
be amplified by merely reducing the pulse width of the long wave component (say 
L). The study on the collision dynamics reveals the interesting behaviour that the 
solitons which split up in the short wave components undergo shape changing collisions 
with intensity redistribution and amplitude-dependent phase shifts. Even though 
similar type of collision is possible in (1+1) dimensional multicomponent integrable 
systems, to our knowledge for the first time we report this kind of collisions in (2+1) 
dimensions. However, solitons which appear in the long wave component exhibit only 
elastic collision though they undergo amplitude-dependent phase shifts. 
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1. Introduction 

One of the main emphasis of current research in the area of integrable systems and their 
apphcations is the study on muhicomponent nonhncar systems admitting sohton type 
solutions [1-17]. In (1+1) dimensions, it has been shown that the muhicomponent bright 
sohtons of the integrable N-coupled nonlinear Schrodinger (CNLS) equations undergo 
fascinating shape changing collisions with intensity redistribution which have no single 
component counterpart [4-8] . This interesting behaviour found applications in nonlinear 
switching devices [12], matter wave switches [13] and more importantly in the context 
of optical computing in bulk media [14, 15]. There is a natural tendency to look for such 
kind of collisions in higher dimensions. From this point of view, we have considered the 
following recently studied integrable coupled (2+1) dimensional ((2+l)D) system, which 
is a two component analogue of the two dimensional long wave-short wave resonance 
interaction (LSRI) system [16], in dimensionless form, 

z{Sif^ + S^^)-S^ + LS^^^O, J = 1,2, (la) 

2 

where the subscripts denote partial derivatives. (Note that in the above two components 
mean two short wave (S) components). The one component {j — 1) version of the above 
equations corresponds to the interaction of a long interfacial wave (L) and a short surface 
wave (S) in a two layer fluid [18]. Also in ref. [19], the existence of dromion like solution 
was established for the j = 1 case. In their very recent interesting work, Ohta, Maruno 
and Oikawa [16] have derived equation (1) as the governing equations for the interaction 
of three nonlinear dispersive waves by applying a reductive perturbation method. Here, 
among these three waves, two waves are propagating in the anomalous dispersion region 
and the third wave is propagating in the normal dispersion regime. In the context of 
long wave-short wave interaction, the first two components can be viewed as the two 
components of the short surface waves while the last component corresponds to the 
long interfacial wave. Note that the presence of the long interfacial wave induces the 
nonlinear interaction between the two short wave components which leads to nontrivial 
collision behaviour as will be shown in this paper. Here onwards wc call equation (1) as 
the 3-wavc LSRI system in which the first two components correspond to short waves 
and the last one is a long wave. 

Apart from deriving the governing equation (1) in ref. [16], Ohta et al have also 
given Wronskian type soliton solutions of a specific type where the components -S''^^^ 
yS'^^) , and L comprise of N solitons, M solitons, and (M + N) solitons, respectively. In 
this context, however, it is of considerable interest to study the collision behaviour if the 
same number of solitons are split up in all the three components and to check whether 
nontrivial shape changing collisions of solitons as in the case of CNLS systems [4, 5] 
occur here also and to look for the possibilities of construction of logic gates based on 
soliton collisions. For the one component case the interaction of two solitons in both 
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short wave and long wave components have been studied in detail in ref. [18] and certain 
interesting features such as fusion and fission processes have been revealed. In this study, 
we consider the multicomponent (2+l)D LSRI system admitting the same number of 
bright solitons in all the three components and obtain the multisoliton solutions. Our 
analysis on their collision properties shows that the solitons appearing in the short wave 
components exhibit a shape changing collision scenario resulting in a redistribution of 
intensity as well as amplitude- dependent phase shift whereas the long wave component 
solitons undergo standard elastic collisions only but with amplitude-dependent phase 
shifts. We also point out that the {N, M,N + M) soliton solutions obtained in ref. [16] 
follow as special cases of the (m, m, m) multisoliton solution obtained here when some of 
the soliton parameters are restricted to very special values. The study is also extended 
to the (n + 1) wave system as well, where n is arbitrary. 

The plan of the paper is as follows: In section 2, we briefly present the bilinearization 
procedure for the three wave system. Multisoliton solution of the three wave system 
is discussed in section 3. Explicit one-soliton and two-soliton solutions are analysed in 
section 4. The asymptotic analysis of the two-soliton solution of the three wave system is 
given in section 5. The interesting collision scenario of two solitons is discussed in detail 
in section 6. Sections 7 and 8 deal with three- and four-soliton solutions, respectively. 
Multicomponent case with j > 2 in equation (1) is studied in section 9. Section 10 is 
allotted for conclusion. 



2. (2-|-l)D bright soliton solutions 

The soliton solutions of equation (1) are obtained by using Hirota's direct method 
[20,21]. By performing the bilinearizing transformations 

where ^''^■'•''s are complex functions while / is a real function, equation (1) can be 
decoupled into the following bilinear equations 

{t{D, + Dy)-Dl){g^^^-f) = Q, J = 1,2, (3a) 

2 

DtD,{f-f) = -2j2i9^^^9^'^*), m 

where * denotes the complex conjugate. The Hirota's bilinear operators D^, Dy and 
are defined as 

DlDlDUa . 6) = ( ^ - ^ ) ( ^ - 77- ) ( ^ - ^ ) aix, y, t)hix', y' , t') 
^ y > \dx dx'J \dy dy'J \dt dt' J ^ ^ ^ ^ ^^=^,^y=y,^t=t,) 

Expanding g^^^^s and / formally as power series expansions in terms of a small arbitrary 
real parameter x, 

9^'^ = X9? + X'9? + ---, J = 1,2, (4a) 

/ =i+xV2+xV4+---, m 
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and solving the resultant set of linear partial differential equations recursively, one can 
obtain the explicit forms of g^'^ and /. Then by substituting their expressions in (2) 
one can write down the soliton solutions. The procedure has been successfully used to 
unearth several interesting properties of soliton collisions associated with CNLS system 
in Refs. [4-8]. We have used a similar procedure here and obtained the one-soliton (1, 
1, 1), two-soliton (2, 2, 2), three-soliton (3, 3, 3) and four-soliton (4, 4, 4) solutions 
explicitly. This can be generalized to the arbitrary m-soliton (m, m, m) solution, in a 
Gram determinant form. Prom this one may claim that in the general case the number 
of sohtons which split up in the short wave components (S"*^^^ and S^'^'') as well as in 
the long wave component (L) are the same. However, we also point out that the (1, 
1, 2), (2, 2, 4) and (N,M,N + M) soliton solutions obtained by Ohta et al [16] can 
be deduced as special cases of our (m, m, m) soliton solution with m = 2, m = 4 and 
m — N + M, respectively, for particular choices of parameters in the solutions. 



3. Arbitrary m-soliton solution 



We first present the general form of (m, m, m) soliton solution for arbitrary m in the 
following Gram determinant form. In order to write down the multisoliton (m-soliton) 
solution of the three wave LSRI system (1), we define the following (1 x m) row matrix 
Csj s = 1,2, (m X 1) column matrices ipj, and 0, j — 1,2, ... ,111, and the (m x m) 
identity matrix /: 



a 



a 



is) 



2 ) ■ ■ ■ ) 



), = (0,0,...,0), 




7 = 



/I 
1 

V 



\ 


1 / 



(5a) 



{5b) 



Here a 



1,2, j = l,2,...,m, are arbitrary complex parameters and rii 



kiX — (ikf + uJi)y + ujit, i = 1,2, . . . ,m, and ki and uji are complex parameters. Then we 
can write down the multisoliton solution of the three wave LSRI system in the form of 
equation (2), with 
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«=1.2, L=-2— log(/). 



(6) 



where 



A 



A I 
-I B 



(7a) 



I <t> 
-I B 0^ 
C, 

in which s denotes the short wave components. Here the matrices A and B are defined 
as 



A 



ki + kj 



, B 



' ' ^, 1,3 = 1,2,. ..,m. {7b) 
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In equation ( l7^ . f represents the transpose conjugate and the real parts of (or 
fcj's) should be chosen as negative quantities in order to obtain nonsingular solutions, 
which are necessary conditions. Sufficiency condition requires the choice of parameters 
such that / is real and nonzero (see below sections 4, 7 and 8 for details in the case of 
m = 1, 2, 3, and 4). 



3.1. Proof of multisoliton solution of the three wave LSRI system 

We now prove that the Gram determinant forms of g^'^^ and / given above indeed satisfy 
the bilinear equations (3). By applying the derivative formula for the determinants, that 



IS, 



dD 

dx 



Edtti j dD 



dx dtti j 

l<jj'<n " 



E 

l<j j'<n 



dx '''' 



la) 



where D 



ail 0.12 

0-21 0-22 



O'ln 
0'2n 



and Aj j is the cofactor of the («,j)^^ element 



and making use of the properties of bordered determinants and also the elementary 
properties of determinants [21,22], the derivatives gx \ fx, ft, fxt, g^z \ gix, fz and fxx, 
where ^ = + ^ j , can be derived as below: 



A I (f) (f)x 

-I B 0^ 



0-10 



fx 



A 



I 
B 0^ 




ft 



s=l 



A I 
-I B 
Cs 



0^ 

^ s 





xt 



E 



A 


/ 




0^ 


-I 


B 


0^ 


-CI 


-0t 























A 


/ 


0- 


EX 




A 


/ 





(f)x 


-I 


B 


0^ 


T 




-I 


B 


0^ 


0^ 





Cs 





+ i 





Cs 














-1 




0t 











/ 





4^xx 




A 


/ 





4>x 




B 


0^ 


0^ 




-I 


B 


0^ 


0^ 




Cs 








+ 





Cs 













-1 







0t 














{8d) 
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A 
-I 



I 4>x 

B {f 




A 
-I 

-4 



I 4> 

B 0^ 




m 



and 





A I 


(t>x 




A 


I 





fxx 


-I B 




+ 


-I 


B 


0^ 




-0t 







-€ 









The conjugate of g^^'' can be written as 

/ 0^ 
B -Cl 




A 
-I 



m 



Substituting for gi^\ g'i\ gi^, /j,, /^.j,, and fz in equation (3a), we find 



A 







/ 

B 
C,, 






0^ 






0a; 

0^ 






A 



I 
B 



A 
-I 



A 
-I 




I 0x 

B 0^ 




A 
-I 




I 

B 



I 

B 



4>x 

0^ 





A 
-I 





0^ 



/ 

5 0^ 




(8z) 



This is nothing but a Jacobian identity and hence g^^^ and / satisfy the first bihnear 
equation (3a). fn a similar way one can also check that the second bilinear equation 
(3b) gives rise to the following Jacobian identity for the Gram determinant forms of g'^^^ 
and /: 



2 

E 

s=l 



A I 

-I B 

-0t 

c. 



0^ 





0^ 

-Cl 






A 
-I 



I 
B 



-E 



A 









/ 
B 



/ 

B 

Cs 



0^ 





0^ 

^ s 



A 



A 

-I 
-0t 



I 
B 0^ 




/ 
B 




0^ 



m 



Thus equations (8i) and (8j) clearly show that the given Gram determinants g^^^ and / 
satisfy the bihnear equations (3), which completes the proof of (6) with (7). 



3.2. {N, M,N + M) soliton solution 

We now point out that the (A^, Af, N + M) soliton solution (for N even) given in ref. [16] 
can be obtained as a special case of the above (m, m, m) soliton solution for the specific 
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choice of parameters af^ — 0,i — 1,2, . . . , N, and — 0,1 — N + 1, N + 2, . . . , m{— 
N + M) along with the parametric restrictions 

a^P - JJ^^ i = l,2,...,N, 



c^f)= J-I^-^^^'^^^l / = 7V + l,iV + 2,...,m(=iV + M). 

In the following sections, we will consider the exphcit cases of m = 1, 2, 3, and 4 soliton 
solutions and the nature of the soliton interactions therein. 



4. One-soliton (1, 1, 1) and two-soliton (2, 2, 2) solutions 



Specializing to the case of m = 1 in equation (6) so that the Gram determinants take 
the form 



9 



All 


1 


e>7i 




-1 












-a? 








An 1 
-1 ^11 



1,2, 



(9) 



where An 



, and Bii — ku 



(1)|2 



\a 



(2)|2^ 



h + K 

explicit one-soliton solution as 



a 



(i). 



~ (log(l + 



J = 1,2, 



One can write down the 

(10a) 
(106) 



where 



r]i = kix — {ikl + u!i)y + Uit, - 
h = kiR + ikii, ui = uiR + iuju- 



(10c) 
(10(i) 



Here a^, af\ Ui and /ci are all complex parameters. In equation (10) the suffixes 
R and / denote the real and imaginary parts, respectively. It may be noted that this 
bright soliton solution is nonsingular only when kiRUOiR < 0, otherwise equation (10) 
becomes singular. In this work, the main focus will be on nonsingular solutions as they 
are of physical importance. The above one-soliton solution can also be rewritten as 



S^^^ = A,VW^e^^"sech (^rj.R + | 
L = -2A;i^sech^ f ?7ii? + — 



1,2, 



(11a) 
(116) 



where 



ViR = kiRX + {2kiRkii - uiR)y uiRt and Aj 



a? 



a 



(1)|2 



-|- a 



(.)|.y 



1,2. 
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The complex quantities Aj\/kiRUJiR, j = 1,2, represent the amphtude of the sohton in 
the S^^^ components whereas thy real quantity 2klj^ gives the amplitude of the soliton 
in the component —L. Note that the complex quantities Ai and A2 satisfy the relation 
|Aip + |y42p = 1, which is a reflection of the fact that the set of equation (1) is 
rotationally symmetric in the {S^^\ S^"^^) space. 

For illustrative purpose, let us obtain the soliton solution for the special choice of 
parameters uJi = —ikf/2. In this case, the above soliton solution (11) becomes 

R 



J = 1,2, 



L = — 2fc^j:jsech^ 



VlR + 



R 



where 



?7i = kix 



E 



-ikfj^ku 



ki = km + ik 



111 



12a) 



[12b) 



:i2c) 



The above soliton solution is nonsingular only when ku < 0, otherwise the parameter R 
in equation f ll2cp becomes complex and the solution (12) becomes singular. Interestingly, 
we observe that by just reducing the width of the soliton in the L component (which 
is proportional to ku) without affecting its amplitude, the soliton in the S^^^ and 
S*^^-* components can be amplified with a proportionate pulse compression, a desirable 
property for a pulse in nonlinear optics. 



4.I. Two-soliton (2, 2, 2) solution 

To obtain the two soliton solution, we take m = 2 in equation (7) and deduce the Gram 
determinant forms as 



Au 


A12 


1 







A21 


A22 





1 




-1 





Bn 


B12 








-1 


B21 


B22 











-a? 


-a? 






An 


A12 
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A21 


A22 





1 


-1 





Bn 


B12 





-1 


B21 


B22 



where A^ 



- — — , and Bij 

ki + k* 



f 



(1) (1)* , (2) (2)* 

a} 'ay + a) 'a} 
{ujj + u*) 



(13) 



then write the explicit form of the (2, 2, 2) soliton solution as 



+ a 



(i). 



,V2 



+ e 



1,2, 



L = 
where 



log(/), 

_j_ QVl+Vl+V2+ri2+R3 



1, 2. We can 

(14a) 
(146) 



'14c) 
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The various quantities found in equation (14) are defined as below: 

■q- = - {iki + uji)y + Uit, i = 1,2, e^^ = ^\ (15a) 

[Ki + Ki) 

{h + k*,y {k, + k*,y {k, + kiy ^'""^ 



{ki - k2 



{k, + kl){k2 + kl) 



e'^^ = 7^^%T7^^^(«2''^'^i2 - ^?^22), J = 1,2, {15d) 



{k2 + kl){ki + kl) 



|2 



R-i 1^1 / \ r \ 

^ = Ti — I — , ,^M, — -T7I^('«ii'«22 - ft:i2/«2i), (15e) 
(fci + kl){k2 + kl)\ki + kl\^ 

( (1) (1)* , (2) (2)*\ 

\a\ 'a\ + a\ ' a\ j 

(^i + ^r) 

The two-sohton solution is characterized by eight arbitrary complex parameters ai' , 
0L\' , 01^2^ , 0L^2^ , ki, k2, uji and 002- The above solution features both singular and 
nonsingular solutions. The nonsingular solution can be obtained by requiring the 
denominator / in (14) to be real and nonzero. The expression (14c) for / can be 
rewritten as 

/ = 2e''^«+^2^ (e(«i+fi2)/2 ^Qgh (r/ifi - 7]2r + (Ri + i?2)/2) + e''^^ cos (r^u - 7]2i + 601) 
+e^^/2 ^Qgi^ ^^^^ ^ ^ ^^/2)) . (15/) 

To get regular solutions, e^^ and e"^^ should be positive which can be obtained only 
for kiROJiR < and k2RUJ2R < 0, respectively. Otherwise, that is for negative values, 
the solution is not regular as in this case -Ri and R2 appearing in the argument of cosh 
in first term become complex. So the condition kjRUjR < 0,j = 1,2 is a necessary 
condition to obtain regular solution. In a similar way, in the third term, the quantity 
R3/2 becomes real and positive for the condition K,ni^22 — I^i2p > 0, as may be seen 
from equation (15e). Still the middle term cos(77i/ — 772/ + Sqi) can lead to a singularity 
as it oscillates between -1 and 1. This can be eliminated by choosing the coefficients 
of the remaining two terms as e('^i"'"^2)/2 _|_ ^R3/2 ^ e''"^, in order to ensure that / will 
not be zero at any point in space and time. The last condition is a sufficient one. As 
an illustration, the interaction of two solitons in system (1) is shown in figure [H The 
parameters are chosen as fci = 1 — 2i, /c2 = 1.5 — 1.05i, ui = —1 — i, uj2 = —1.3 — 0.5i, 
a^i^ = 2, a'^^ = ap"* = 1, aj^"* = 0.01. One observes that the solitons in the S^'^^ and 
S*^^-* components undergo shape changing (energy redistribution) collisions while there 
is only elastic collision in the L component. More details are given in section 6 below. 

4-2. (1, 1, 2) soliton solution of Ohta et al 

Now we show that the (1, 1, 2) soliton solution obtained by Ohta et al [16] is a special 
case of the above two-soliton (2, 2, 2) solution (14). Specifically, for the special choice 



Bright solitons in multicomponent LSRI system 



10 



of the parameters 0:2^^ = — 0, the above two-sohton solution becomes 

1 



L = -2 



(log(/)) , 



where 



J _ ]^ _|_ ^Vi+rii+Ri _|_ gr?2+T?2+-f^2 _|_ ^m+rit+m+ri^+Rs 

The various other parameters defined in equations (14) now take the forms 



Jill 



(fcl+fci*)' 

„(2) 



«22 



(16a) 
(16&) 
(16c) 

(16(i) 

(16e) 

(16i) 
(16^) 

(a;i+a;i*)' '^^^ " (a;2+a;2*)- ^^^''^ 
Solution (16a-h) is nothing but the (1, 1, 2) soliton solution obtained by Ohta tt al in 
ref. [16] when the parameters in (16a-h) are further restricted to the special choice 



(^2 + ^1)' 



\k\ - A;2pKllK22 



-Q;[^^fi;22(fc2 - /^i) 
(A:2 + A;|)(A;i + A;|)' 



= 0, 



I (1)12 
Fi T 



/«22 



;i) ^ (/Ci + fct)(fci + fe^) 

^ (A;2 - k^) 



and a 



(2) ^ (fe2 + fc2*)(fc2 + fei*) 

(A;2 - -fci) 



(16z) 



5. Asymptotic analysis of the two soliton solution (14) of the three wave 
system 



We now consider the collision properties associated with the general two-soliton solution 
(14) of the three wave system. For this purpose we carry out the analysis, for kj^i > 0, 



1,2. Also we choose > 



and 



kiRkii 



> 



1^2 R 
^IR 



for convenience. 



Similar analysis can be performed for other choices of kjjis and cOjRS also by keeping 
kjji > 0, UjR < 0, which is the necessary condition for nonsingular solutions. We 
now define the soliton wave variables as rjiR = kmx + {2kiRkii — uJiR)y + uJiRt and 
ViR = k2RX + {2k2Rk2i — oj2R)y + 0J2Rt- In the limit x,y ^ ±00 and a fixed t the 
two-soliton solution (14) takes the following asymptotic forms, 
a) Before collision (limit y — > —00): 
(i) Soliton 1 (rjiR ~ 0, 772R ^ -00): 

/ -5(1) \ f Al-\ 



\J kiRUJiR sech {r]iR + ^] e*"" 



-2klji sech^ 



Ri 



(17a) 



(176) 



Bright solitons in multicomponent LSRI system 
where 



( A\-\ 



((A;i + A;t)(cui + cut))V2- 



(ii) Soliton 2 {r]2R ~ 0, r]iR oo): 
/ 5(1) \ I Al-\ 



\ '''' J 



\Al- J 



where 



\J k2R^2R sech ^r72R + ^ 
(i?3 — Ri] 



-(Ri+Ri)/2 



\ '''' ) 



The various quantities in the above expressions are defined in equation (15). 

b) After coUision (hmit y — > oo): 
(i) Sohton 1 {rjiR ~ 0,772ii ^ oo): 



V kiRUJiR sech ^r/iR + 



(-R3 — -R2; 



where 



Lc.-2kl^ sech^ 



{i?2+-R3)/2 



((A;i + A;i*)(a;i+cUi*))V2- 



(ii) Sohton 2 (r72i? - 0, r^iR -00): 



\ '''' J 



■2RUJ2R sech {r]2R + ^] e*'"' 



where 



(R 
V2R+-^ 

\Al- J 



-R2/2 



{{k2 + kl){uJ2+Ujl)yl^' 
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Note that in all the above expressions + |^2^|^ = 1; J = 1)2. Our above 

analysis reveals the fact that due to collision the amplitude of the colliding solitons, 
say si and $2 in the S^^^ and S**^^^ components, change from (^A\~ ,A\~^ VkiRtoiR and 
{Al'',Al') ^/k2RUJ2R to ^/kiRUiR and v^W^, respectively. 

Here the superscripts in Aj^'s with z,j = 1,2 denote the solitons si and S2, while 
the subscripts represent the components S^^^ and S''-^-* and the " ± " signs stand for 
"x, y —>■ ±cxd" . In addition to this change in the amplitudes, the solitons also undergo 
amplitude- dependent phase shifts due to the collision and they can be determined 
straightforwardly from the above asymptotic expressions. From equations ((17) and 
(19)) and equations ((18) and (20)), one can easily check that the phase shift suffered 
by the soliton si (say $i) = — Phase shift of soliton S2 (say —$2 = '^'i) = ^ and is 
given by 

^^ (R^-R^-R^) ^ (21) 

where Ri, R2 and -R3 are as defined in equation (15) and depends on the amplitudes. 



6. Soliton Interaction 



Now it is of further interest to analyze the interaction properties of the solitons depicted 
in figure [T] for the specific set of values of the parameters given in section 4.1. Figure 
[1] shows typical spatial collision of two solitons for t = — 4 corresponding to the exact 
expression (14). The interesting collision scenario depicted in figure [1] clearly indicates 
that there is a redistribution of intensity among the two S^^^ components resulting in 
an enhancement (suppression) of intensity of solitons S2 (si) in the S^^^ component and 
a suppression (enhancement) of soliton S2 (si) in the S**^^-* component. The solitons 
also undergo amplitude-dependent phase shifts along with this energy redistribution. 
However, the solitons appearing in the long wave component (L) exhibit the standard 
elastic collision as shown in the third figure of figure [T] though the phase shift here 
is also amplitude-dependent. Interestingly, if the parameters are so chosen such that 
the condition -nj = -fay is satisfied, there occurs only elastic collision in all the three 

components S^^^ and L. The underlying collision dynamics can be well understood by 
using the asymptotic analysis of the two-soliton solution (14) discussed in section 5 and 
is further described below. 



6.1. Collision behaviour of solitons in the short wave components 

The asymptotic analysis presented in the previous section also results in the following 
expressions relating the intensities of solitons si and S2 in the S^^^ and S**^^^ components 
before and after interaction, 

= |t;| Vrr, ^,J = 1,2, (22a) 

where the superscripts j± represent the solitons designated as Si and S2 at "x, y ±00" . 
The expression for the transition intensities for the solitons in the short wave components 



Bright solitons in multicomponent LSRI system 



13 



can be written down using the results in equations fllTd . (118 d . fllQcl) and (120 cp as 



1|2 



IT' 



Ai 



'2 1 2 



A1A2I 
A1A2I 



|1 - 
^21 



1,2, 



A2 



(225) 
(22 c) 
(22(i) 



In general |Tj| 7^ 1 and so an intensity (energy) redistribution of the solitons in the 
S^^^ and S**^^-* components occurs as shown in figure [TJ One can notice that the standard 
elastic collision takes place for the specific parametric choice 



it 
.(1) 



(2) 

(2), dS 



and hence = l^^"*^!^, i,j = 1, 2, for this choice. However the two colliding solitons 

si and S2 suffer amplitude-dependent phase shifts $1 and $2, respectively, as given in 
equation fl^ . 



6.2. Collision scenario in the long wave component 

In the L component, there occurs only elastic collision for any parametric choice. This 
is evident from the asymptotic analysis, vide equations (17b), (18b), (19b) and (20b). 
One finds that the amplitudes of the solitons Si and $2 before and after interaction are 
the same which are —2k'l^ and —2^27?, respectively, while there occurs an amplitude- 
dependent phase shift as given by equation fl2T|) . 



6. 3. Shape changing collisions and Linear fractional transformations 

It is instructive to notice that the intensity redistribution in the short wave components 
characterized by the transition matrices (equation (22)) can also be viewed as a linear 
fractional transformation (LFT). To realize this, we re-express the amplitude changes 
in the short wave components of soliton Si after interaction as 

A\+ = VCiiA\- + VCi2A\-, (23a) 
TC2iA\- + TC22A\-. (236) 



Here 



where 



A\ 



a 



[al 0.2 + oi\ 0.2 )\P-2 "2 + ^2 '^2 



(2)^(2)*w^(l)^(l)* 



„(2)^,(2)*^ 



C12 
C21 
C22 



("2 "2 )(^1 - ^2) + (^2 "2 )(^1 



(a2^^a?^*)(^^2 + ^^2), 
(tt?^a2^^*)(^^2 + ^^2), 

\afaf*){u^ + ^) + {afaf*){yj^ - ^2) 



-1/2 



(23 c) 
(23(i) 

(23e) 
(23/) 
(23^) 
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ih - + kl){co2 + cu*)(aS')a^')* + af^af *)] . (23z) 

Note that the coefficients Cj/s, i,j = 1,2, are independent of a^^^s and so of and 
Al~ , that is the a parameters of sohton Si. From equations (123 ap and (123 6p . 



where p\~2 = 4t^5 which the superscripts represent the underlying sohton and 
the subscripts represent the corresponding short wave components. Thus the state 
of Si before and after interaction is characterized by the complex quantities pl~^ and 
pW, respectively. The direct consequence of the above LFT representation is the 
identification of a binary logic using soliton collisions as in the case of CNLS equations 
[6, 14, 15] and hence the LFT can be profitably used to construct logic gates associated 
with the binary logic. A similar analysis can be made for the soliton S2 also. 

7. Three-soliton (3, 3, 3) solution 

From the general form (7), and restricting m = 3, one can write down the explicit 
three-soliton (3, 3, 3) solution as 



a 



f 

g'72+'72+'?3+'54j _|_ QV3+V3+m+i5j _|_ ^VS+Va+m+^Sj _|_ g»?^+'72+»73+'57j _|_ g''?l+»?2+'73+% 

+ 

+ 



/ 



where 



Here 



/ 

^V2+V2+V3+ri3+Vl +T3j 

+ J , J = 1,2, (24a) 



J — ]_ _l_ gVi+iit+Ri _|_ g'?2+'72+^?2 _|_ ^V3+vi+R3 _|_ Qm+V2+^u) _|_ Qvl+m+^to 

_|_ Qm+V3+S2Q _|_ QVi+m+S20 -)- g'?2+r?3+(S30 _|_ g'?2+'?3+'530 _|_ g'*?! +'?^ +»?2+r?2 +-^4 

_|_ gf?l+'7t+'?3+'7|+-R5 _|_ gf?2+»75+r?3+'?|+-fi'.6 _|_ Qm+-ni+V2+V3+'ri0 _j_ gm+'73'+'?3+f?2+'^10 

_j_ ^V2+'ri2+Vl+V3+'^20 _j_ g''?2+»?2+'?l+'?3+T-2*o _|_ g»?3+'7|+f?l+''y5+f"30 _|_ ^m+Vs+Vt+m+r.^o 

_j_ QVl+Vt+V2+V2+V3+V3+R7 (246) 

rji = kiX — {ik^ + u!i)y + ujit,i = 1,2,3, (24c) 
Jii _ ~ k2){a[^^K2i - a^2^Kn) _ {ki - k3){a^^^K3i - a^^^Kn) 



{k, + kl){kl + k2) ' {ki + kl){kl + k 



3) 



5„. {ki-k2){a[^^K22-ai^^ni2) 5,, {k2 - k3){a^2^ K32 - a^3^ 1^22) 



gJ3j _ ^ 1 ^z, Z -L^/ g04j 



(A;i + fc2*)(^2 + A;2*) ' {k2 + k*){k*2 + k^) 



■em 
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"mm 



(^2 - A:3)(q;2'^^^33 - a?K23) 
{kl + k2){kl + k3) 

{kl - k3){a'f^K32 - 0(3^ K12) 
(ki + k*2)(k*2 + k3) 



(fci-A-,)(np\v23- 
{ki + kl){k2 + kl) 

{k2 - ki){k3 ~ ki){k3 - k2){kl - k*) 
{kl + ki){kl + k2){kl + k3){k*2 + ki){kl + A;2)(A;2* + ^3) 

X a^i\K2lK32 - K22«3l) + a2^H'^12«31 " «32fi;il) + OLi\KiiK22 " fi;i2fi;2l) 
(fc2-fcl)(fc3-fci)(fc3-fc2)(fc;-fcn 

{kl + ki){kl + k2){kl + k3){kl + A;i)(A;3* + A;2)(A;3* + /C3) 

X Cxf\K33n2l - K31K23) + 02^^31 ft^lS " ^n/tss) + Q!3-'\/«23/«ll - /«13/«2l) 

{k2-ki){k3-ki){k3-k2){kl-kl) 

{k*2 + A;i)(A;2* + A;2)(A;2* + A;3)(A;3* + A;i)(A;| + k2){k*3 + k3) 

Qi?^(/«22«33 - «23«32) + Qi2^^(«13«32 " «33«12) + Oi3\'^l2l^23 " «22/«13) 



kn 



m 



1,2,3, 



«12 



,520 _ 



ki + kl 



D<530 



A;i + A;^ ' 



^23 
/C2 + /Sq 



(24(/) 
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The explicit form of L can be obtained by substituting the expression for / in 
L = -2 — (log/). Here a'-^\ aj^^ aj^^ af\ af\ ki,k2,h, tui, wa, and 

are the twelve complex parameters which characterize the above three-soliton solution. 
Following the arguments of ref. [8] (see equations (28) and (29) there) and the discussion 
in section 4 one can show that the necessary conditions for nonsingular solution are 

which is automatically taken care by the choice kj^uojn < 0,j = 1,2,3. One can also 
easily show that the inequality 

g(/?l+ii6)/2^ g(iJ2+i?5)/2^ g(i?3+iJ4)/2^ g(^?7)/2 ^ 4inaX (^e^^WR+T-iOR) ^ g{S20R+r20R) ^ ^iSsOR+rum)^ (^24:h) 

is the sufficient condition in order to ensure that the solution is regular. Typical shape 
changing collision of the three-soliton solution is shown in figure [2] for the parametric 
choices = 0.2 + 0.3i, fcs = 0.6 + 0.4i, fcg = 0.7+0.2z, tui = -0.5 + OAi, 002 = -0.7+O.li, 
UJ3 = -0.3 + 0.3i, af ^ = 0.5 - i, a^^^ = 0.5 + i, a^^^ = 0.3 + 0.2i, af ^ = 0.39 + 0.2i, and 
a'^^ = a^^^ = 1. The above three soliton solution (24) represents the interaction of three 
solitons and their collision scenario can be well understood by making an asymptotic 
analysis following the procedure given in section 5 for the two soliton solution. 

We have identified from the asymptotic analysis that for the three interacting 
solitons (say, si, S2 and S3), as in the case of CNLS equations [6,8], the total transition 
amplitude of a particular soliton (say Si) can be expressed as the product of two 
transition amplitudes which result respectively during the first collision of si with S2 
and during the collision of the outcoming soliton (say s'l) with soliton S3. In a similar 
manner the net phase shift acquired by a particular soliton (say si) during the complete 
collision process is equal to the addition of phase shifts experienced by that soliton 
during its cascaded collisions with S2 and S3, respectively. Thus the analysis clearly 
shows that the multi-soliton collision process indeed occurs in a pair-wise manner in 
the multicomponent (2-|-l)D LSRI system and there exists no multiparticle effects. The 
details are similar to the CNLS system [6,8] and so we do not present them here. 



8. Four-soliton (4, 4, 4) solution 

Again to obtain the explicit four-soliton solution, we substitute m = 4 in the Gram 
determinant form (7) and obtain an expression involving exponentials. Since it is too 
lengthy, we do not present the explicit form here. However we note that the four-soliton 
solution is characterized by sixteen complex parameters, kj, uj, ctf'\ ctf \ j = 1,2,3,4. 
The nonsingular solution results for the choice kjRUjR < 0, j = 1,2,3,4. One can 
check that kjRUJjR < 0, j = 1,2,3,4, are the necessary conditions for the existence of 
nonsingular solution and the sufficient condition can be obtained following the procedure 
mentioned in sections 4 and 7. 
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1. (2, 2, 4) soliton solution of Ohta et al 

the above discussed four-soliton solution, we make the choice a"^'' = a^^^ = af^ 

(2) 



— 0, and also introduce the parametric restrictions 

(1) ^ ih + ki){h + k;){h + ki){h + ki) 

{k,-h){k,-k,){k,-k,) ' ^'^"^ 

(1) ^ {k2 + kl){k2 + kl){k2 + kl){k2 + kl) 

{kl-k2){k,-k2){k,~k,) ' ^^'"^ 

(2) ^ {h + kl){k^ + ki){k^ + ki){k^ + ki) 

3 ^k,-k,){k2-k,){k,-k,) ' ^ ' 

(2) _ (^4 + kl){k^ + A:*)(fc4 + fc*)(fc4 + fc*) 

(fc3-A:4)(A:2-A:4)(fci-fc4) ' ^ ' 

then one can show that it is exactly equivalent to the (2, 2, 4) soliton expression given 
by Ohta et al [16]. This can be verified by expanding the determinant form with 
the above parametric restrictions and comparing it with the expanded version of the 
(2, 2, 4) solution of ref. [16]. The interaction of solitons for the above special case of the 
four-soliton solution is shown in figure [3] for the choice of parameters ki = 0.5 — 0.2i, 
k2 = 0.4 + O.li, ks = 0.3 - OAi, k^ = 0.4 + 0.6i, ui = -0.5 + 0.4i, 002 = -0.7 + O.li, 
UJ3 = —0.3 + 0.3i, uji = —0.2 + 0.2i. This is similar to the interaction shown in ref. [16]. 

9. Soliton solutions of {n + l)-wave system 

We now extend our study to obtain multisoliton solutions of the multicomponent system 
with arbitrary {n + 1) waves, in which we consider n short wave components and single 
long wave component. The {n + l)-wave system in this case is given by 

z(5p^ + 4^'))-5(i) + L5(^') = 0, j = l,2,...,n, (26a) 

n 

Lt = 2j2\S^'^l- (266) 
(i) One soliton solution: 



Following the procedure discussed in section 4, we can obtain the one soliton solution 

as 

^''= l+',n..,t.n ^ J = l,2,...n, (27a) 
L= -2^(log(l + e^^+''^*+«)), (276) 



where 



r^, = kix - {tkl + uji)y + uit, = ^''^^ ' ^ \ (27c) 



fii) Two soliton solution: 



Similar procedure results in the two soliton solution for the multicomponent case with 
arbitrary {n + 1) waves whose expression can be obtained from equations (14) by just 
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allowing j to run from 1,2, ... ,n and redefining ku as ku — — X]j=i('^l''^Q;P^*)/('^«+'^r)' 
i,l — 1,2. Now it is straightforward to extend the bilinearization procedure of obtaining 
one and two soliton solutions to obtain multisoliton solutions as in the ID integrable 
CNLS equations [8]. Similarly, three- and four soliton solutions of equation (26) can be 
obtained by suitably redefining n^s with i,l = 1,2,3 and i,l = 1,2,3,4, respectively, 
and fixing the upper limit of the index j, corresponding to the short wave components 
as n. The multisoliton solution of the multicomponent case (26) can be written down 
from equation (6) by allowing s to run from 1 to n and redefining the column matrix 
ipj as 

i^j = (af\af\...,af''^^ . 

Likewise the proof of the multisoliton solution of the multicomponent system also follows 
the three wave system discussed in section 3. 

10. Conclusion 

To conclude, we have obtained explicitly the multi bright plane soliton solutions of 
recently reported physically interesting integrable (2+1) dimensional (n+l)-wave system 
by applying Hirota's bilinearization procedure. We have also presented the results in 
a Gram determinant form for the multisoliton solutions of the multicomponent LSRI 
system along with the necessary proof. We observe that the solitons in the short wave 
components can be amplified by merely reducing the pulse width of the long wave 
component. The study on collision dynamics shows that the solitons appearing in the 
short wave components undergo shape changing collisions with intensity redistribution 
and amplitude-dependent phase shift. This gives the exciting possibility of soliton 
collision based computing in higher dimensional integrable systems also. However, the 
solitons in the long wave component always undergo elastic collision. 
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Figure 1. Shape changing coUision of soUtons in the three wave system with the 

(1) (2) 

parametric choice ^^jj ^ -jri- Here the two soliton solution (14) is plotted for a fixed 

value of t and the associated parameters are given in the text (below equation (15)). 
Note that intensity redistribution occurs only in the 5'*^^^ and 5*^^^ components, while 
elastic collision only occurs in the L component. 
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Figure 2. Shape changing colhsion of the three sohtons of the three wave system, 
equation (24). The chosen soHton parameters are given in the text. Again note that 
intensity redistribution occurs only in the S'^^^ and S'*^^-' components, while elastic 
collision only occurs in the L component. 




Figure 3. A special case of the four soliton collision with two solitons in S^^^ and S^'^^ 
components and four solitons in L component, equation (25). Note that this collision 
scenario is similar to the (2, 2, 4) soliton interaction depicted in figure 2 of ref. [16]. 



